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Abstract 

This paper constructs continuously self-similar solution of a spherically sym- 
metric gravitational collapse of a scalar field in n dimensions. The qualitative 
behavior of these solutions is explained, and closed-form answers are provided 
where possible. Equivalence of scalar field couplings is used to show a way to 
generalize minimally coupled scalar field solutions to the model with general 
coupling. 
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I. INTRODUCTION 



Clioptuik's discovery of critical phenomena in gravitational collapse of a scalar field |I| 
sparkled a surge of interest in gravitational collapse just at the threshold of black hole 
formation. Discovery of critical behavior in several other matter models quickly followed 
0-^. While perhaps the presence of critical behavior in gravitational collapse is not in itself 
surprising, some of its features are, in particular the conclusion that black holes of arbitrary 
small mass can be formed in the process. Moreover, the critical solution often displays 
additional peculiar symmetry - the so-called self-similarity, and serves as an intermediate 
attractor for near-critical solutions. 

Study of critical phenomena also throws new light on the cosmic censorship conjecture. 
The formation of strong curvature singularity in critical collapse from regular initial data 
offers a new counterexample to the cosmic censorship conjecture. 

Much work has been done, and general features of critical behavior are now understood. 
However, there is a distinct and uncomfortable lack of analytical solutions. Due to the 
obvious difficulties in obtaining solutions of Einstein equations in closed form, most of the 
work seems to be done numerically. One of the few known closed form solutions related to 
critical phenomena is Roberts solution, originally constructed as a counterexample to the 
cosmic censorship conjecture [Q, and later rediscovered in the context of critical gravitational 
collapse P|T(]|]. It is a continuously self-similar solution of a spherically symmetric gravita- 
tional collapse of a minimally coupled massless scalar field in four dimensional spacetime. 



For review of the role self-similarity plays in general relativity see ||TT 



This paper searches for continuously self-similar, spherically symmetric scalar field solu- 
tions in n dimensions. They might be relevant in context of the superstring theory, which 
is often said to be the next "theory of everything" , as well as for understanding how critical 
behavior depends on the dimensionality of the spacetime. Roberts solution would be a par- 
ticular case of the solutions discussed here. These solutions provide reasonably simple toy 
models of critical collapse, although they are not attractors |]12[. Some qualitative properties 
of the self-similar critical collapse of scalar field in higher dimensions have been discussed in 
| 13| . Here we aim at finding explicit closed-form solutions. 

Second part of this paper deals with extension of minimally coupled scalar field solutions 
to a wider class of couplings. It is shown how different couplings of scalar field are equivalent, 
and several particular models are examined in details. Procedure discussed here can be 
applied to any solution of Einstein-scalar field equations. 



II. REDUCED ACTION AND FIELD EQUATIONS 



Evolution of the minimally coupled scalar field in n dimensions is described by action 
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(1) 



plus surface terms. Field equations are obtained by varying this action with respect to field 
variables g^i, and 0. However, if one is only interested in spherically symmetric solutions 
(as we are), it is much simpler to work with reduced action and field equations, where this 
symmetry of the spacetime is factored out. 
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Spherically symmetric spacetime is described by the metric 



ds^ = d'j^ + r^dto'^ 



(2) 



where c?7^ is the metric on a two-manifold, and du'^ is the metric of a (n — 2)-dimensional 
sphere. Essentially, the spherical symmetry reduces the number of dimensions to two, with 
spacetime fully described by the two-metric dj'^ and two-scalar r. It can be shown that the 
reduced action describing field dynamics in spherical symmetry is 



'S'sph OC 



-7 d^x R[-f] + {n-2){n- 3)r-^ ( Vr 



2(V0)^ 



(3) 



where curvature and differential operators are calculated using two-dimensional metric 7as; 
capital Latin indices run through {1,2}, and stroke | denotes covariant derivative with 
respect to two-metric 'Jab- Varying reduced action with respect to field variables jab, 
and 0, we obtain Einstein-scalar field equations in spherical symmetry. After some algebraic 
manipulation they can be written as 



Rab -{n~ 2)r-V| 



AB 



(4a) 



n-3) (Vr)2 - 1 



□0 + (n - 2)r~S^'^0,Ar 



B 



(4b) 
(4c) 



As usual with scalar field, equation is redundant. 



III. iV-DIMENSIONAL GENERALIZATION OF ROBERTS SOLUTION 

We are interested in generalization of Roberts solution to n dimensions. To find it, we 
write the metric in double null coordinates 

dj^ = -2e-^''^''^dudv, r = -up{z), (p = (p{z). (5) 

The dependence of metric coefficients and (p on z = —v/u only reflects the fact that we 
are looking for continuously self-similar solution, with z being scale-invariant variable. We 
turn on the influx of the scalar field at the advanced time w = 0, so that the spacetime is 
Minkowskian to the past of this surface, and the initial conditions are specified by continuity. 
Signs are chosen so that 2; > 0, p > in the sector of interest (m < 0, f > 0). With this 
choice of metric, Einstein-scalar equations become 

(n - 2) [p"z + 2a' p - 2a' p' z] = -2pz0'^ (6a) 
2p{a"z + a') + (n - 2)p"z = -2p#'^ (6b) 
{n - 2) [p" - 2a' p'] = -2p0'^ (6c) 
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(n-3) 



+ p"pz = 0, 



(6d) 



'pz + (n - 2)(f)'p'z -^{n- ^)<\> p = 0. 



(6e) 



Prime denotes the derivative with respect to z. Combining equations (|6a| ) and we 
obtain that a = const. By appropriate rescahng of coordinates, we can put a = 0. Then 



(n - 2)p" = -2p0 



/2 



(7a) 



(n-3) 



P Z-PP+- 



+ p"pz = 0, 



(7b) 



^ + (n-2)^-i(n-4)^-i = 0. 



(7c) 



For further derivation we will assume that n > 3, as the case n = 3 is trivial. Equation (^ 
can be immediately integrated 



(8) 



Substituting this result back into equation ([7aD, we get the equation for p only 



^•"O n-4 



n - 2 



(9) 



It is easy to show that equation (D is equivalent to equation ([751) . No surprises here, since 
the system (^) was redundant. Combining both equations we get first integral of motion 



P z- PP+- 



z 



n— 3 
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{n-2){n-3)' 



(10) 



which contains only first derivatives of p, and for this reason is simpler to solve than either 
one of equations (|7E| ) and (^). Equation (|10|) is generalized homogeneous one, and can be 
solved by substitution 



X = - Inz, 
2 
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p = y/zy{x), p' = ^'"^{y + y) 



(11) 



where dot denotes the derivative with respect to new variable x. With this substitution, 
equations (p!OD and (||) become 



= y'^ -2 + ciy 



-2{n-3) 



(12) 



2coy 



-(n-2) 



(13) 
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where we redefined constant 



= ^ > 0. (14) 



The above equation (KT) for y formally describes motion of a particle with zero energy in 
potential 

V{y) = 2-y'-c,y-''-''-'\ (15) 



so we can tell qualitative behavior of y without actually solving equation ([1^). 

Initial conditions are specified by continuous matching of the solution to Minkowskian 
spacetime on surface v = 0. Since on that surface r ^ 0, the value of y = r j ^/—uv starts 
from infinity at x = — oo, and rolls towards zero. What happens next depends on the shape 
of the potential. If there is region with V{y) > 0, as in Fig. |I|, y will reach a turning point 
and will go back to infinity as x = oo. If V{y) < everywhere, as in Fig. there is nothing 
to stop y from reaching zero, at which point singularity is formed. Finally, if V{y) has 
second-order zero, as in Fig. ^, y will take forever reaching it. 

Of course, variables separate, and equation ([T^) can be integrated 



x = ±l , '^^ +C2. (16) 

ly^-2 + ciy-2("-3) 



Plus or minus sign in front of the integral depends on the sign of the derivative of y. Initial 
conditions imply that initially y comes from infinity towards zero, i.e. its derivative is neg- 
ative, and so we must pick the branch of the solution which started out with a minus sign. 
Constant C2 corresponds to a coordinate freedom in choice of origin of x, while constant Ci 
is a real parameter of the solution. 

Unfortunately, the integral cannot be evaluated in a closed form for arbitrary n. But if 
the integral is evaluated, and we can invert it to get y as a function of x, the solution for r 
is obtained by using definitions (|TT]) and (^. The solution for is obtained by integrating 
relations (ITSf) or f|8|). 



IV. CRITICAL BEHAVIOR 

One-parameter family of self-similar scalar field solutions in n dimensions we constructed 
above exhibits critical behavior as the parameter ci is tuned, much as Roberts family does 
in four dimensions. In this section we investigate black hole formation in the collapse. 

In spherical symmetry, existence and position of the apparent horizon are given by van- 
ishing of (Vr)^ = 0, which translates to p' = 0, or ?/ + ?/ = in our notation. Therefore, at 
the apparent horizon we have 

y'-y' = ciy-^(--3) -2 = 0, (17) 

and the black hole is formed if value of y reaches 
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!/i„ =(!)". (18) 

As we have discussed above, depending on the value of ci, values of field y either reach 
turning point and return to infinity, or go all the way to zero. The critical solution separates 
the two cases, and is characterized by potential Viy) having second order zero, i.e. Viy^) = 
V'{y^:) = at some point Differentiating expression ( |T5| ) for potential V{y), we see that 
it has second order zero at 

(19) 

if and only if the value of the constant Ci is 
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n — 3 



2.^ 



n-2 



(20) 



If the value of parameter ci is less than critical, Ci < c^, the value of y turns around at the 
turning point, and never reaches point of the apparent horizon formation, which is located 
in the forbidden zone, as illustrated in Fig. This case is subcritical evolution of the field. 
If Cl > cl, the value of y reaches point where apparent horizon is formed, and proceeds to 
go to zero, at which place there is a singularity inside the black hole. This supercritical 
evolution is illustrated in Fig. ||. 

The mass of the black hole formed in the supercritical collapse is 

1 1 

M = -rAH = -^u^/z^yAH■ (21) 

It grows infinitely if we wait long enough, and will absorb all the field influx coming from 
past infinity. This happens because the solution is self-similar, and creates a problem for 
discussing mass scaling in the near-critical collapse. Cut and glue schemes avoiding 
infinite black hole mass are temporary means to lift this problem. However, the real answer 
to determining if the critical solution is intermediate attractor and calculating mass-scaling 
exponent is to make a perturbative analysis of critical solution. Similarity to Roberts solu- 
tion suggests that the results for four dimensional spacetime [0 can be applied to higher 
dimensions as well. 



V. PARTICULAR CASES 



In this section we consider several particular cases for which the general solution (|T6D 
is simplified. Particularly important is n = 4 case, which is the already familiar Roberts 
solution. 



A. n = 3 

As we have already mentioned, for n = 3 the only self-similar scalar field solution of the 
form (^ is trivial. To see it, note that equation ( [75| ) implies that p" = if n = 3, and so 
p = az + P and r = av — (3u. From equation ( [7a]) it then follows that = const. The 
spacetime is fiat. 
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B. n = 4 



Integration (|^) can be carried out explicitly 

dy 



X = — 



C2 



2 



- 1 + VZ/^ - + ci 



+ C2, 



(22) 



and the result inverted 



(23) 



to give the solution in the closed form 



P 



'e2'=2 1 - ci 

— + Z + 



2e2^2 



(24) 



By appropriately rescaling coordinates, we can put e^'^^ = 2. After redefining parameter of 
the solution p = (ci — l)/4, the solution takes on the following simple form 



p = \Jl + z — pz'^, r = \Ju'^ — uv — pv"^. 



(25) 



Scalar field (f) is reconstructed from equation (P) 



_2 1 v^r+4^ 

6 = cop - 



2 1 + 2;— pz'^ ' 



(26) 



to give 



, 2pz-l 
arctanh „ + const 



2 



v/TT4^ 
2pz - 1 + VI + 4p' 

2pz-l - y/1 + 4j9 



+ const. 



(27) 



The critical parameter value is p* = 0, and for p > the black hole is formed. The critical 
solution is 



Vu'^ — 



1 



In 



1 - 



(28) 



C. n = 5,6 



Integral (16) can be written in terms of elliptic functions, which becomes apparent with 
the change of variable y = y~^ 



X 



dy 



yy/1 -2y + ciy"^-'^' 
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(29) 



and the solution y{x) is given implicitly. However, integrals corresponding to critical solu- 
tions simplify, and can be taken in terms of elementary functions for n = 5,6. It happens 
because potential factors since it has second order zero, therefore reducing the power of y 
in the radical by two. The results of integration for critical solutions are 



n = o : X 



y^dy 



{y' - 1) 



- arcsmh h/^.U ^ arctanh ( ^MH] + ' arctanh ( ) (30) 




n = 6 : X = — 



y^dy 



{y' - 1) ^y' + y' + l 

arcsinh ( V2(?/^ + o) ) + TTTS arctanh ( ^ = | . (31) 



2 V 27 2v^ \V2,/^r^ + 



The dependence y(x) is still given implicitly. The critical value of the parameter is 32/27 
for n = 5, and 27/16 for n = 6. 

D. Higher dimensions 



For higher dimensions, integral (|T^) can not be taken in terms of elementary functions, 
so one has to be content with the solution in the integral-implicit form, or do numerical 
calculations. 



VI. GENERAL SCALAR FIELD COUPLING 

In this section we discuss in detail how solutions of minimally coupled scalar field model 
can be generalized to much wider class of couplings. The fact that essentially all couplings 
of a free scalar field to its kinetic term and scalar curvature are equivalent has been used in 
the past |TB|,|TB[ to study scalar field models with non-minimal coupling. In particular, this 



idea has been applied to extend four- dimensional Roberts solution to conformal coupling 
ITl and Brans-Dicke theory . 



A. Equivalence of Couplings 

Suppose that the action describing evolution of the scalar field in n-dimensional spacetime 

is 

3 = 77- [V^d\ \F{(f))R - (!>)'] (32) 

plus surface terms, where the couplings F and G are smooth functions of field 0. Also sup- 
pose that the signs of couplings F and G correspond to the case of gravitational attraction. 
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We will demonstrate that this action reduces to the minimally coupled one by redefinition 
of the fields g^^ and 0. First, let's introduce new metric g^j^y that is related to old one by 
conformal transformation 



— 
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g = n\ 



-9, 



(33) 



and denote quantities and operators calculated using g^y by hat. Scalar curvatures calculated 
using old and new metrics are related 



R = Q^R + 2{n - l)QnQ - n{n - 

field gradients 

(V0)' = fi'(V0)^ 
Writing action (^) in terms of metric gfj_^, we obtain 
1 



(34) 



(35) 



S 



gd\ F {n' R + 2{n - l)nnn - n{n - l){Vny } -Gn\V(f)y . (36) 



By choosing conformal factor to be 



n-2 



(37) 



the factor in front of the curvature R can be set to one. Substitution of definition of Q into 
the above action, and integration by parts of □ operator yields 



S 



gcTx 



G n-1 F' 



n-2 F"^ 



^ (V0)^ 



(38) 



The kinetic term in action (^Sj) can be brought into minimal form by introduction of a new 
scalar field 0, related to the old one by 



2(V0)' 



F^^2F^]^^^^- 



(39) 



Thus, we have shown that with field redefinitions 



1 



'G F'^\ 



g^iu F ^9lJ.i^: 



V2J \F n-2 F^ J 
the generally coupled scalar field action ( |32D becomes minimally coupled 
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(40) 
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gd\ R-2{V(py 



(41) 



This equivalence allows one to construct solutions of the model with general coupling ( |32D 
from the solutions of minimally coupled scalar field by means of inverse of relation (^OD, 
provided said inverse is well-defined. However, there may be some restrictions on the range 



9 



of SO that field redefinitions give real and positive-definite metric g^,y. This means 
that not all the branches of the solution in general coupling may be covered by translating 
minimally coupled solution. Technically speaking, the correspondence between solutions of 
minimally coupled theory and generally coupled theory is one-to-one where defined, but not 
onto. 

However, one has to be careful making claims about global structure and critical behavior 
of the generalized solutions based solely on the properties of the minimally-coupled solution. 
The scalar field solutions encountered in critical collapse often lead to singular conformal 
transformations, which could, in principle, change the structure of spacetime. 



B. Examples 



To illustrate the above discussion, we consider two often used scalar field couplings as 
examples. They are non-minimal coupling and dilaton gravity. 



1. Conformal coupling 



Non-minimally coupled scalar field in n dimensions is described by action 
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gd''x (1 - 2^0")i? - 2(V0)^ 



(42) 



where ^ is the coupling parameter. Coupling factors are F = 1 — 2,^0^, G = 2 and so the 
field redefinition (H0|) looks like 



1 - 2^02 



^ - arcsin 



arcsinh 



(43) 



where 



1 n-2 



4 n - 1 

Particularly interesting is the case of conformal coupling ^ = because field redefinition 

1 



(44) 



arctanh 



(45) 



can be inverted explicitly to give the recipe for obtaining conformally coupled solutions from 
minimally coupled ones. It is 



tanh 



(46) 
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cosh^/("-2) 



1 - 2^,02)2/("-2) 

In particular, four dimensional Roberts solution becomes 



-\/3 tanh 



1 2pz-l 
—pz arctann , 

x/3 



(47) 



(48) 



ds = cosh 



1 2pz-l 
—pz arctann , 

^/3 VI +4p 



in conformally coupled model. This last expression was considered in [|T^]. 



(49) 



2. Dilaton gravity 

Another useful example is dilaton gravity described by the action 

1 



S 



167r 



-gd^e-'f i? + 4(V0) 



(50) 



Substituting coupling factors F = e G = — 4e '^'^ into relationship (^OD, one can see that 
the scalar field redefinition is a simple scaling 



n-2 



ln-2 



and metrics differ by exponential factor only 

Qf,^ = exp 



n-2 



In particular, four dimensional Roberts solution becomes 

6 = arctanh f ^ , 

vrT4^' 



ds"^ = e^*^ {-2du dv + (m^ - uv - pv^) doo^} 



in dilaton gravity. 



(51) 



(52) 



(53) 



(54) 
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VII. CONCLUSION 



We have searched for and found continuously self-similar spherically symmetric solutions 
of minimally coupled scalar field collapse in n dimensional spacetime. For spacetime dimen- 
sions higher than three they form one-parameter family and display critical behavior much 
like the Roberts solution. Qualitative picture of field evolution is easy to visualize in analogue 
with particle traveling in potential of upside-down U shape. Unfortunately, the solutions in 
dimensions higher than four can only be obtained in implicit form. Critical solutions are 
in general simpler than other members of the family due to the potential factoring. Strong 
similarity between Roberts solution and its higher dimensional generalizations allows one to 
conjecture that these higher dimensional critical solutions are not attractors either. Absence 
of non-trivial self-similar solution in three dimensions raises question whether scalar field 
collapsing in three dimensions exhibits critical behavior at all. Perhaps further numerical 
simulations will answer it. 

We also use equivalence of scalar field couplings to generalize solutions of minimally 
coupled scalar field to much wider class of couplings. For often-used cases of conformal 
coupling and dilaton gravity the results are remarkably simple. Some results of , applied 
for single scalar field only, become trivial in view of this coupling equivalence. 

However, the question of critical behavior of these generalized solutions is complicated by 
the fact that the conformal factor f2 relating metrics (|33|) for minimally and generally coupled 
solutions may be singular. In the simplest case of non-singular conformal transformation 
(i.e. when coupling F is bounded and lower bound is greater than zero) global properties 
of the minimally coupled solution are preserved, and all important features of near-critical 
collapse carry over on the generalized solution unchanged. If conformal transformation is 
suspected to be singular, more careful study of global properties of generalized solution (^Of ) 
is necessary. 
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FIG. 1. Subcritical field evolution. 
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FIG. 2. Supercritical field evolution. 
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FIG. 3. Critical field evolution. 
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